Master equation for a kinetic model of trading market and its analytic solution 
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We analyze an ideal gas like model of a trading market with quenched random saving factors for 
its agents and show that the steady state income (m) distribution P(m) in the model has a power 
law tail with Pareto index v exactly equal to unity, confirming the earlier numerical studies on this 
model. The analysis starts with the development of a master equation for the time development of 
W)' P{m). Precise solutions are then obtained in some special cases. 
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' The distribution of wealth among individuals in an economy has been an important area of research in economics, for 
O more than a hundred years. Pareto 1] first quantified the high-end of the income distribution in a society and found 
' it to follow a power-law P{ra) mr''^^"\ where P gives the normalized number of people with income m, and the 

a I exponent i^, called the Pareto index, was found to have a value between 1 and 3. 
Considerable investigations with real data during the last ten years revealed that the tail of the income distribution 
' indeed follows the above mentioned behavior and the value of the Pareto index v is generally seen to vary between 1 
, and 2.5 j^l^Q. It is also known that typically less than 10% of the population in any country possesses about 40% 
Q of the total wealth of that country and they follow the above law. The rest of the low income population, in fact the 
O majority (90% or more), follow a different distribution which is debated to be either Gibbs j^, ||5| or log- normal 0|. 

Much work has been done recently on models of markets, where economic (trading) activity is analogous to some 
\ scattering process [1, IE S S • We put our attention to models where introducing a saving factor for the 

^ • agents, a wealth distribution similar to that in the real economy can be obtained |E ■ Savings do play an important 
' role in determining the nature of the wealth distribution in an economy and this has already been observed in some 
■ recent investigations |12| |. Two variants of the model have been of recent interest; namely, where the agents have the 
' same fixed saving factor and where the agents have a quenched random distribution of saving factors . While 
the former has been understood to a certain extent (see e.g, 0,^3)' ^'^^ argued to resemble a gamma distribution 
[T^ . attempts to analyze the latter model are still incomplete (see however, Further numerical studies 0| of 

\ time correlations in the model seem to indicate even more intriguing features of the model. In this paper, we intend 
"Y^ ' to analyze the second market model with randomly distributed saving factor, using a master equation type approach 
' similar to kinetic models of condensed matter. 



■ II. THE MODEL 
C " 

Q , The market consists of N (fixed) agents, each having money mi{t) at time t {i = 1,2, .. . ,N). The total money 
M (= iTiilt)) in the market is also fixed. Each agent i has a saving factor (0 < Ai < 1) such that in any 

• .—I , trading (considered as a scattering) the agent saves a fraction Xi of its money rrii (t) at that time and offers the rest 

^ ' (1 — \i)mi{t) for random trading. We assume each trading to be a two-body (scattering) process. The evolution of 

^ \ money in such a trading can be written as: 



m,(t + 1) = Kmi{t) + [(1 - \i)m,{t) + {I - \j)mj{t)] , (1) 
mj{t + 1) = Xjmj{t) + (1 - ey) [(1 - \i)m,{t) + (1 - Xj)mj{t)] (2) 
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where each rrii > and etj is a random fraction (0 < e < 1). Typical numerical results for the steady state money 
distribution in such a model is shown in Fig. 1(a) for uniform distribution of (0 < < 1) among the agents. 



III. DYNAMICS OF MONEY EXCHANGE 



We will now investigate the steady state distribution of money resulting from the above two equations representing 
the trading and money dynamics. We will now solve the dynamics of money distribution in two limits. In one case, 
we study the evolution of the mutual money difference among the agents and look for a self-consistent equation for 
its steady state distribution. In the other case, we develop a master equation for the money distribution function. 



A. Distribution of money difference 

Clearly in the process as considered above, the total money (mi + mj) of the pair of agents i and j remains constant 



while the difference Arriij evolves as 



(Amy)t+i = (mj-mj)t+i = ^ ^' + ^ {Amij)t+(^ '^' 2 ) + + (3) 

Numerically, as shown in Fig. 1, we observe that the steady state money distribution in the market becomes a power 
law, following such tradings when the saving factor A^ of the agents remain constant over time but varies from agent 
to agent widely. As shown in the numerical simulation results for P{m) in Fig. 1(b), the law, as well as the exponent, 
remains unchanged even when =1/2 for every trading. This can be justified by the earlier numerical observation 
0, for fixed A market (Ai = A for all i) that in the steady state, criticality occurs as A ^ 1 where of course the 
dynamics becomes extremely slow. In other words, after the steady state is realized, the third term in |j2Jl becomes 
unimportant for the critical behavior. We therefore concentrate on this case, where the above evolution equation for 
Amy can be written in a more simplified form as 

(Amy )t4.i = Qfy (Amy )t + ^^(mi + mj)t, (4) 

where aij = + Aj) and [3ij = — Aj). As such, < a < 1 and — ^ < /3 < i. 

The steady state probability distribution D for the modulus A = |Ato| of the mutual money difference between 
any two agents in the market can be obtained from Q in the following way provided A is very much larger than the 
average money per agent = M/N . This is because, using eqn. Q), large A can appear at t+1, say, from 'scattering' 
from any situation at t for which the right hand side of eqn. 10} is large. The possibilities are (at t) nii large (rare) 
and rrij not large, where the right hand side of eqn. becomes ^ (a^ + /?y)(Ay)t; or rrij large (rare) and not 
large (making the right hand side of eqn. Q becomes ^ (a^ — /3ij)(Ay)t); or when rm and rrij are both large, which 
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is a much rarer situation than the first two and hence is neghgible. Then if, say, is large and rrij is not, the right 
hand side of I^J becomes ~ (a^ + Pij){Aij)t and so on. Consequently for large A the distribution D satisfies 

D{A) = J dA' D{A') {S{A - (a + /3)A') + 5{A - {a ~ /3)A')) 



where we have used the symmetry of the (3 distribution and the relation aij + (3ij = Ai, and have suppressed labels 
J, J. Here (...) denote average over A distribution in the market. Taking now a uniform random distribution of the 
saving factor A, p{\) = 1 for < A < 1, and assuming D{A) ^ A~^^+'^^ for large A, we get 

l = 2y" dAA^ = 2(l + 7)-\ (6) 

giving 7=1. No other value fits the above equation. This also indicates that the money distribution P{m) in the 
market also follows a similar power law variation, P{m) ~ m~^^~^'^'> and v = ^. We will now show in a more rigorous 
way that indeed the only stable solution corresponds to = 1, as observed numerically El- 



B. Master equation and its analysis 

We now proceed to develop a Boltzmann-like master equation for the time development of P(rn,t), the probability 
distribution of money in the market. We again consider the case e^j = ^ in J^l and Q and rewrite them as 

where 



""-{i^)^ .-4(1±A). (8) 

Collecting the contributions from terms scattering in and subtracting those scattering out, we can write the master 
equation for P{m,t) as (cf. llj) 



P{m,t + At) - P{m,t) = {J drui J dnij P{mi,t)P{Tnj,t) 

X {[(5({^ m}i — m) + S{{A m}j — m)] — [5{fni — m) + S{mj — m)]}} 
= ildm.ldm,Pim.,t)Pim„t) 

X [S{iifmi + nij — m) + S{fi^ nii + j-ij-nij — m) — 5{mi — m) + 5[mj — m)]). (9) 

The above equation can be rewritten as 

dP{m,t) 
di 

which in the steady state gives 



P{m,t) — { / drrii / drrij P{mi,t)P{mj,t) 6{fj,]'mi + iJ,j nij — m)), (10) 



P{m) = ( / drrii j drrij P{mi)P{mj) 5{iifmi + jij nij — m)). (11) 



Writing rriiii'l = xm, we can decompose the range [0, 1] of x into three regions: [0, k], [k, 1 — k'] and [1 — k', 1]. 
Collecting the relevant terms in the three regions, we can rewrite the equation for P{m) above as 



, , m , „ / xm\ „ / m(l — x) \ , 

P(m) = (^- / dxPi — )pi ^ _ M ) 

- {^\P(-)^ rdyPiy)+p(^)i^ f^dyPiy)+ f'dxP 



(12) 
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where the result appHes for k, and k! sufBciently small. If we take tti^I/k, m^l/re' and k, k' ^ (m — > oo), then 
Assuming now as before, P(ra) — A/m}^'^ for m ^ oo, we get 

1 = {{pi+r + {^l-r) ^ j j d^,+di,-p{i,+)qif,-) [i^i+r + , (w) 

as the ratio of the third term in to the other terms vanishes like (ttik)"'^, [mn')^'^ in this limit and p{fJ-^) and 
are the distributions of the variables /i^ and fi^ , which vary uniformly in the ranges 1] and [0, ^] respectively 
(cf. eqn ©). The i,j indices, for /i+ and /i~ are again suppressed here in 114() and we utilise the fact that /i^ and 
/i~ are independent for i ^ j. An alternative way of deriving Eqn. H14(l from Eqn. is to consider the dominant 

terms (oc x^^ for r > 0, or oc ln(l/a;) for r = 0) in the x — + limit of the integral m'^'^'^'^^ P{m) exp{—mx)dm (see 
Appendix A). We therefore get from Eqn. p4|l . after integrations, 1 = 2/(i^ + 1), giving = 1. 



IV. SUMMARY AND DISCUSSIONS 

In our models 0, S 13 ' consider the ideal-gas- like trading markets where each agent is identified with a gas 
molecule and each trading as an elastic or money-conserving (two-body) collision. Unlike in a gas, we introduce a 
saving factor A for each agents. Our model, without savings (A = 0), obviously yield a Gibbs law for the steady-state 
money distribution. Our numerical results for various widely distributed (quenched) saving factor A showed 01 H 
that the steady state income distribution P{m) in the market has a power-law tail P{m) ~ m~'-^+''-' for large income 
limit, where v ~ 1.0. This observation has been confirmed in several later numerical studies as well [Tsl Il6||. Since 
Q{m) = P[m)dm can be identified with the inverse rank, our observation in the model with v = 1 suggests that 
the rank of any agent goes inversely with his/her income/wealth, fitting very well with the Zipf's original observation 
[TtI I. It has been noted from these numerical simulation studies that the large income gro up people usually have 
larger saving factors |3- This, in fact, compares well with observations in real markets The time correlations 

induced by the random saving factor also has an interesting power-law behavior J 6,] . A master equation for P(m, i), 
as in lO, for the original case (eqns. ^ and Q) was first formulated for fixed A (Ai same for all i), in 13] and solved 
numerically. Later, a generalized master equation for the same, where A is distributed, was formulated and solved in 

Us.]. 

We have formulated here a Boltzmann-type master equation for the distributed saving factor case in Q and jSJ. 
Based on the observation that even in the case with e = 1/2 (with A distributed in the range < Ai < 1, Ai Aj), in 
and the steady state money distribution has the same power-law behavior as in the general case and shows 
the same Pareto index, we solve the master equation for this special case. We show that the analytic results clearly 
support the power-law for P[m) with the exponent value v — 1. Although our analysis of the solution of the master 
equation is for a special case and it cannot be readily extended to explore the wide universality of the Pareto exponent 
as observed in the numerical simulations of the various versions of our model 0, ITsj , let alone the quasi- universality 
for other v values as observed in the real markets ^] , the demonstration here that the master equation admits 
of a Pareto- like power law solution (for large m) with v = 1, should be significant. 

Apart from the intriguing observation that Gibbs (1901) and Pareto (1897) distributions fall in the same category 
of models and can appear naturally in the century-old and well-established kinetic theory of gas, our study indicates 
the appearance of self-organized criticality in the simplest (gas-like) models so far, when the stability effect of savings 
is incorporated. This remarkable effect can be analyzed in terms of master equations developed here and can also be 
studied analytically in the special limits considered. 
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APPENDIX A: ALTERNATIVE SOLUTION OF THE STEADY STATE MASTER EQUATION ([TTll 



Let Sr{x) = dmP{m)m''+'' exp{-mx)- r>0,x>0. If P(m) = A/m^+'', then 



5*^(2;) = A dm ^exp(— ma;) 
"'0 



A if r > 

Aln ( - j if r = 0. 



From eqn. we can write 

/•oo 

5,. (a;) ( 

Jo Jo 

/>oo 

dmi Amp"^(exp(-mi/i+a;) (/^+) } 



' dmi I dmj P{mi)P{mj){mi^f +mjjij )'^^'' exp[— (mi/x^ + mjiij )x]) 
JO 



dmj P{mj){exp(—mjfij x)) 



+ I dmj Am^j ^(exp(-mj^^ x) (^^ Y^^) 



dmi P{mi){exp{~mifif x)) 



or, 



(Al) 



(A2) 



Sr{x) = J^^ d^il p{^f) (^J^ dm^Aml ^ exp(-mj^+x)^ 

+ dfiJq{^J)(^J dmj Amy^ exp{-mj^j,Jx)^ {fijY^'' , 
since for small x, the terms in the square brackets in ljA2|l approach unity. We can therefore rewrite HA3|I as 



Sr{x) ^ 2 



dii+{^i+y+'-Srixfi+)+ / dfi-i^i-y+'-Srixfi-) 

'0 



(A3) 



(A4) 



Using now the forms of Sr{x) as in (|A1|) . and collecting terms of order x ^ (for r > 0) or of order ln(l/.T) (for 
r = 0) from both sides of (|X4|) . we get ((nil . 
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